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Abstract 

Using the canonical Arnowitt-Deser-Misner Hamiltonian formalism, a “first law of mechanics” is 
established for binary systems of point masses moving along generic stable bound (eccentric) orbits. 
This relationship is checked to hold within the post-Newtonian approximation to general relativity, 
up to third (3PN) order. Several applications are discussed, including the use of gravitational self¬ 
force results to inform post-Newtonian theory and the effective one-body model for eccentric-orbit 
compact binaries. 
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I. INTRODUCTION 


The direct observation of gravitational waves would have a tremendous impact on physics, 
astrophysics and cosmology [1], Binary systems composed of compact objects, namely black 
holes and neutrons stars, are highly promising sources of gravitational waves [2], The orbital 
dynamics of these systems can be investigated by means of various approximation methods in 
general relativity, such as post-Newtonian (PN) theory [3-6], the perturbative gravitational 
self-force (GSF) approach [7-10] and the effective one-body (EOB) framework [11-13]. These 
approximation methods are complementary to fully nonlinear numerical simulations of the 
late inspiral and merger of compact binaries. Recent years have seen a spur of activity at the 
multiple interfaces of these various techniques. Such cross-cultural studies have contributed 
to improving our knowledge of the two-body dynamics and wave emission [14]. 

Some of this progress stems from the “Erst law of binary mechanics” of Le Tiec et al. [15] 
(hereafter Paper I), a variational formula that relates small changes in the bodies’ masses to 
those of the total mass-energy and angular momentum of a binary system moving along a 
circular orbit. This variational relationship is a particular case, valid when one assumes the 
existence of a global helical Killing held and a point-particle model for the compact objects,^ 
of the generalized hrst law of mechanics of Friedman et al. [16]. More recently, Blanche! et 
al. [17] have extended to spinning point particles the first law of Paper I, for spins aligned 
or anti-aligned with the orbital angular momentum. 

The hrst law of Paper I was successfully applied to (i) the determination of the numerical 
values of some previously unknown PN coefhcients (at 4PN, 5PN and 6PN orders) that enter 
the expressions for the binding energy and angular momentum as functions of the circular- 
orbit frequency [15], (ii) the derivation of the exact expressions for the binding energy and 
orbital angular momentum at linear order in the mass ratio [18], and (iii) the computation 
of the frequency shift of the Schwarzschild innermost stable circular orbit (ISCO) frequency 
induced by the conservative piece of the GSF [18, 19]. Moreover, it was used in Refs. [19, 20] 
to compute one of the potentials that enters the EOB ehective metric, exactly, at linear order 
in the mass ratio. 

The hrst laws of Refs. [15-17] were established for compact binaries moving along circular 
orbits. In this paper, we shall extend the hrst law of Paper I to binary systems of nonspinning 
compact objects moving along generic stable bound (eccentric) orbits. There are multiple 
reasons to do so. First, although most stellar-mass compact binaries would have completely 
circularized by the time they enter the observable frequency band of ground-based detectors 
such as Advanced LIGO, Advanced Virgo and KAGRA, there are scenarii where eccentricity 
ehects could become observable and would give access to much interesting physics [21-26]. 
Second, eccentric inspirals with extreme mass ratios are promising sources for a future mHz- 
band gravitational-wave antenna in space such as the proposed eLISA mission [27-31]. Third, 
eccentric orbits give access to new degrees of freedom in the EOB model [32]. 

The remainder of this paper is organized as follows. The hrst law of binary mechanics for 
eccentric-orbit binaries is derived in Sec. II, while some of its consequences are explored in 
Sec. III. This relationship is checked to hold true, in the context of the PN approximation, 
up to 3PN order, in Sec. IV. Finally, Sec. V is devoted to applications, such as the use of 
GSF results to inform PN theory and the EOB model for eccentric-orbit compact binaries. 
Throughout this paper we use “geometrized units” where G = c = 1. 

^ Mathematically, the approximation of an exactly closed circular orbit translates into the existence of a 
helical Killing vector field, along the orbits of which the spacetime geometry is invariant. 
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II. DERIVATIONS OF THE FIRST LAW 


In this section, building upon the canonical Arnowitt-Deser-Misner (ADM) Hamiltonian 
formalism (Sec. II A), we derive a hrst law of mechanics for binary systems of compact objects 
moving along eccentric orbits. The hrst proof relies on an orbital averaging (Sec. IIB), while 
the second proof makes use of action-angle type variables (Sec. IIC). 


A. Canonical Hamiltonian of point-particle binaries 

Our starting point is the application of the ADM canonical formulation of general relativ¬ 
ity [33] to a binary system of nonspinning compact objects, modelled as point particles with 
constant masses rua (with a = 1, 2), as recently reviewed in Refs [4, 34]. The orbital dynamics 
of this binary system is assumed to derive from an autonomous Hamiltonian Pa; 'ma), 

where Xa(t) and Pa{t) are the canonical positions and momenta of the two particles, dehned 
in some hxed gauge. Such a Hamiltonian has been computed, up to 4PN order, ^ using the 
ADM transverse-traceless (TT) gauge [35-37]. The phase-space coordinates Xa(t) and Pa(t) 
obey Hamilton’s equations 

_ OH ^ OH 

Tl ) Pa 
^Pa 

where Euclidean three-vectors are denoted in bold font, and an overdot stands for the deriva¬ 
tive d/dt with respect to coordinate time t. 

In general relativity, asymptotically hat spacetimes possess ten conserved quantities, given 
as surface integrals at spatial inhnity, that are associated— via Noether’s theorem—with the 
continuous symmetries of the Poincare group [38, 39]. These quantities are the total energy 
H, total linear momentum P, total angular momentum L, and boost vector K = G — P t, 
where G is the center-of-mass vector. These act as “generators” of time translations, spatial 
translations, spatial rotations, and Lorentz boosts, respectively. For binary point-particle 
spacetimes, the global Poincare symmetry is realized by the Poincare algebra satished by 
the ten generators 77, P, L and K, viewed as functions on the two-body phase space (x^, Pa) 
[40]. If the coordinate system used to dehne the canonical Hamiltonian manifestly respects 
the Euclidean group, such as the usual ADM-TT gauge conditions, then 77(xa, p^; ma) must 
be translationally and rotationally invariant, which implies [40] 

P = ^ Pa , L = Xa X Pa 
a a 

Thereafter, we shall consider only the relative motion of the binary system with respect 
to the center-of-mass frame, which is dehned by the condition G = 0, implying P = G = 0. 
This restriction decreases by six the number of degrees of freedom of the dynamical system. 
Because of translational invariance, the center-of-mass Hamiltonian 77(r, p; ma) is a function 
of the relative position r = xi — X 2 and the relative momentum p = pi = —p 2 only. Thus, in 
the center-of-mass frame, the conserved orbital angular momentum simply reads L = r x p, 
such that the orbital motion is conhned to the coordinate plane orthogonal to L = L L and 
spanned by r{t) and p{t) at any given time t. Introducing polar coordinates in this plane. 


^ Starting at 4PN order, the binary Hamiltonian involves some nonlocal-in-time contribution, such that H 
becomes a functional of the phase-space variables Xa(<) and Pa{t). 


( 2 . 2 ) 
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such that r = (r cos y9, r sin y9, 0) with r = |xi — X 2 I the separation and ip the orbital phase, 
we have L = = const. Therefore, Hamilton’s equations (2.1) reduce to 


dH 


r = 


Pr = 


dH 


if = 


dH 

'dL 


L = 


dH 


= 0 


(2.3) 


dpr dr dL dip 

and the center-of-mass Hamiltonian H{r,pr, L] rria) is independent of the cyclic cordinate ip. 

Henceforth, we shall consider only bound (and stable) orbits, for which the radial motion is 
bounded: r G [rmin, ?"max], where rmin and r^ax denote the coordinate separation at periastron 
and apastron, respectively. These two turning points of the radial motion correspond to the 
two largest (real, positive, and finite) roots of the equation pr = —dH{r,pr, L;ma)/dr = 0, 
evaluated at pr = 0. Since the binary dynamics is conservative, the radial motion of a bound 
orbit must be periodic in time, such that r{t + P) = r{t) and pr{t + P) = Pr{t) for all time 
t, where P denotes the coordinate time period of the radial motion. 


B. First law from an orbital averaging 


Under infinitesimal changes hr, 6pr, SL and of the phase-space coordinates and mass 
parameters, the Hamiltonian undergoes a variation 


6H 


dH ^ dH 


£ £ T 


dH 

drua 


5ma 


(2.4) 


where to simplify the notation we do not indicate that the partial derivatives with respect 
to r, Pr, L and rua are computed while keeping the other variables fixed. 

For two point particles interacting only through (Einsteinian) gravitation, it was shown 
in Ref. [17] that the partial derivatives of the canonical ADM Hamiltonian with respect to 
the particle’s masses, holding the phase-space coordinates fixed, are simply given by 


dH _ 
drUa 


(2.5) 


where Tait) is the proper time elapsed along the worldline of particle a. The functions Zaif) 
are refered to as the redshift variables, e.g. [41-50]. The proof of Eq. (2.5) involves three main 
steps: (i) the construction of a Fokker Lagrangian [51] for the binary system by eliminating 
the gravitational field degrees of freedom in the total (matter plus held) Lagrangian, (ii) a 
perturbative redehnition of the particle’s coordinate positions to remove all accelerations and 
higher-order time derivatives from this Fokker Lagrangian [52], (iii) a Legendre transform 
to obtain an ordinary Hamiltonian; see Sec. HI of Ref. [17] for more details. These steps are 
analogous to those followed to compute, in the context of the post-Newtonian aproximation, 
an ordinary (Fokker-type) Hamiltonian for two point masses.^ 

If the changes Sr{t), Spr{t) and 6L map two neighboring solutions of the binary dynamics, 
then Hamilton’s equations (2.3) must be satished, and we obtain 


6M = —pr 6r + r 6pr + ip 6L + Za dnia , 

a 


( 2 . 6 ) 


^ After imposing the usual ADM-TT gauge conditions, the field variables hj^- and hj^- are eliminated from 
a Routh functional i?(xa, Pa, , hj^-) constructed from the total (matter plus field) ADM Hamiltonian, 
yielding a (matter only) higher-order Hamiltonian iJ(x^, , Xq , , • • •) that can be reduced to an ordinary 

Hamiltonian through a canonical transformation; see e.g. Refs. [35, 53, 54]. 
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where we used Eq. (2.5) and the fact that, “on shell,” the Hamiltonian is numerically equal 
to the ADM mass M. Since r(t), (fit), Pr(t), Za(t), 6r{t) and 6pr(t) are not constant, the 
relation (2.6) is not particularly useful. 

However, we may get rid of these orbital variations by taking the average of Eq. (2.6) over 
one radial period of the motion. Introducing the notation (/) = A /(t) dt for the time 
average of any function / (recall that P is the coordinate time period of the radial motion), 
and using the fact that M, L and iria are constants of the motion, we get 

5M = (0) 5L + (r 5pr — Pr 5r) + ^^{za) ^rua . (2.7) 

a 

This variational relationship is clearly reminiscent of the first law of mechanics established 
in Paper I for circular motion [see Eq. (1.1) there], with the constant circular-orbit frequency 
and the constant redshifts replaced by their orbit-averaged counterparts (0) and {za). These 
averaged quantities carry clear physical interpretations. Indeed, 

1 /“f $ 

I fP 1 rPa p 

i^a) =pj^ ^ P / ’ 

where $ = 27r -I- A$ is the accumulated azimuthal angle (or orbital phase) per radial period, 
with A<h the relativistic periastron advance, and Ta is the proper time period of the radial 
motion of particle a. Notice that P also coincides with the proper time period of the radial 
motion for a distant inertial observer, where the geometry is essentially flat. 

The radial contribution in Eq. (2.7), which vanishes for circular motion, can also be given 
a simple physical interpretation. Integrating by parts and using the periodicity of Pr(t) Sr(t), 
with period P, we obtain 

(r Spr — Pr Sr) = (r 6pr + Pr Sr) = {S{rpr)) = n SR , (2.9) 

where n = 2% /P is the radial frequency and R = S Pr dr = ^ r'’max radial action 

integral. (The last step in Eq. (2.9) made use of the change of variable f —)■ r in the integral.) 
Finally, introducing the notation oo = (p) for the average azimuthal frequency, we obtain 
the variational relationship 



(2.8a) 

(2.8b) 


SM = u SL + nSR + ''^^{za) SrUa ■ ( 2 - 10 ) 

a 

This is our first law of mechanics for binary systems of nonspinning compact objects moving 
along nearby eccentric orbits. Equation (2.10) generalizes to generic (bound) orbits the first 
law of Paper I, previously established for circular motion. In order to emphasize the structure 
(time period or angular period) x (variation of the “conjugate variable”), Eq. (2.10) may 
also be written as 

PSM = ^SL + 2T^SR + ^TaSma. ( 2 . 11 ) 

a 

The occurrence of the action integrals L = ^ § p^dp and R = ^ f Prdr suggests an alter¬ 
native proof of the first law of binary mechanics, based on the use of action-angle variables. 
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C. First law from action-angle variables 


In the previous section, we studied a 4-diniensional dynamical system with a Hamiltonian 
ip, Pr, p^', rria) that depends on two external parameters rria. In this section, we shall now 
consider an extended, 8-dimensional phase space with 4 generalized coordinates (r, p, Ti, T 2 ) 
and 4 conjugate momenta {Pr,Pip,PTi,PT 2 )- particular, the proper times elapsed along 
the worldlines of the particles now play the role of two additional (non-compact) coordinates. 
Therefore, in addition to Eqs. (2.3) we have the extra canonical equations 

dH . dH 

'^a 7 ^ ) pTa “q 

dPr, dTa 

Since the Hamiltonian does not depend explicitly on the proper times Ta, the momenta pr^ 
must be constants of the motion. Comparing (2.12) with (2.5), it is clear that Pr^ = ma¬ 
llow, this 8-dimensional dynamical system possesses 4 hrst integrals = (M, L, mi, m 2 ) 
that are independent (for non-degenerate orbits) and in involution, i.e., which have vanishing 
Poisson brackets.^ Therefore, this dynamical system is completely integrable [55]. For bound 
orbits, the motion in phase space is bounded in the r and p directions; however it is not 
bounded in the Ta directions, such that the Liouville-Arnol’d theorem does not guarantee the 
existence of action-angle type variables [55]. Nevertheless, a generalization of this theorem 
by Fiorani et al. [56] ensures the existence of generalized action-angle variables {qa, Ja) for 
completely integrable dynamical systems with non-compact level sets; see Secs. H A and B 
of Ref. [57] for a nice summary, given in the language of symplectic geometry. 

Thus, for our dynamical system we can dehne the generalized action variables [56, 57] 


( 2 , 12 ) 




J'.n — 




y Pr dr = R, 

jp^dp = L, 

1 


(2.13a) 

(2.13b) 

(2.13c) 


Since the Hamiltonian H{r,pr, L, rua) = M does not depend on the coordinates p, Ti and T 2 , 
the radial momentum pr can be expressed solely as a function of r and the hrst integrals Pa- 
Hence the action variables Ja = {Jr, J<f, Jn, Jt 2 ) are functions of the hrst integrals only. The 
generalization of the Liouville-Arnol’d theorem for non-compact level sets [56] guarantees 
that these relationships can be inverted, yielding the expressions PaiJg) of the hrst integrals 
in terms of the action variables. 

Then, the complete solution of the Hamilton-Jacobi equation for the action S{t, Sa, Pa) = 
—M t + W{sa; Pa) of the two-body system can easily be found by separation of the variables 
Sa = {r,p,Ti,T 2 ), yielding the following expression for Hamilton’s characteristic function: 


W{Sa;Pa) 


/ r 

Pr{r'-, Pa) dr' + Lp + ^ 


m„Ta 


(2.14) 


From the characteristic function W we dehne a generating function G{sa, Ja) = W (s^; Pa{,Jg)) 
that generates a Type H canonical transformation, yielding a new Hamiltonian H' expressed 


^ This is obvious from the triviality of the first integrals Pa = (i?,p,^,Pri,Pr 2 )- 
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in generalized action-angle variables (g^, Ja)- Since G does not depend explicitly on time t, 
we have Ja) = Hisa,Pa)- The generalized angle variables can be computed from the 

relation qa = dG/dJa- On the other hand, pa = dG/dsa is already satished from Eq. (2.14). 
In generalized action-angle variables, Hamilton’s equations take on the simple form 

OH' ■ dH' 

Qa o 7 — ^a 5 Ja ^ 0 

JJa JQa 

Because the Hamiltonian H'{Ja) does not depend on the generalized angles ga(t), the actions 
Ja are constants of the motion, as expected from the fact that they are functions of the first 
integrals Pa only. Thus the fundamental frequencies Ua are also constants of the motion, such 
that qait) = qa,o +(^at, where qa,o are four constants that reflect a choice of “initial phase.” 
The wide class of coordinate transformations that leave the oJa unchanged is discussed, e.g., 
in Ref. [57] (see also [58, 59]). While Ur and are true angular frequencies, the quantities 
cJt-i and are dimensionless. 

Now, by varying the Hamiltonian H' and using the equations of motion (2.15), we imme¬ 
diately get bH' = 'Yha^o^ Then, using Eqs. (2.13) and the fact that H' = M “on shell,” 
we obtain the variational relation 

6M = ojr SR + oj^bL + biJig . (2.16) 

a 

Moreover, since the fundamental frequencies associated with the (generalized) angle variables 
qr, q,p and g^-a read Ur = 27 i/P = n, = ^/P = u and ojr^ = Ta/P = (za), we recover the 
hrst law (2.10). Notice also that the averaging over one radial period is already taken care 
of when using action-angle type variables. 

Finally, we discuss the special case of circular orbits, which are characterize by a vanishing 
radial action integral: R = 0. For such degenerate orbits, the constants of the motion M and 
L are no longer independent, and the motion becomes simply periodic with constant angular 
frequency uj. In that limit, an interesting coordinate-invariant relationship is given by the 
reduced periastron advance K = uj/n = $/(27r) as a function of the circular-orbit frequency 
oj (see, e.g., Refs. [54, 59-65]), computed from the zero-eccentricity limit oi K = —OR/dL, 
evaluated at fixed M, mi and m 2 . 


(2.15) 


III. CONSEQUENCES OF THE EIRST LAW 

We now explore some straighforward consequences of the first law for eccentric-orbit bina¬ 
ries. In particular, in section HI A we derive various partial differential equations relating the 
coordinate-invariant relationships M{u),n,ma), L{oj,n,ma), R{u),n,ma) and {za){uj,n,ma), 
and we establish a “first integral” associated with the variational law (2.10) in section HIB. 


A. Partial differential equations 


For given particles’ masses ma, the action integrals L and R uniquely specify an orbit up 
to initial conditions. Using (L, R, mi, m 2 ) as independent variables, the variational first law 
(2.10) thus yields the following set of partial differential equations: 


dM _ dM _ dM 
dL ^ ’ dR ^ ’ dma 


(Za) ■ 


(3.1) 
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In words, the partial derivatives of the ADM mass M with respect to the generalized action 
variables L, R, mi and m 2 are simply given by the fnndamental freqnencies u, n, (zi) and 
(^ 2 )- However, it is often more convenient to parameterize the motion using the two frequen¬ 
cies u and n instead of the action integrals L and R. Using (cn, n, mi, m 2 ) as independent 
variables, Eq. (2.10) is equivalent to the relationships 


dM 


duj 

dM 


dn 

dM 


dma 


dL dR 


= u 


= UJ 


= u 


doj 

dL 

dn 

dL 

dma 


n 


n 


duj ’ 
dR 

dn ’ 
dR 

^ dma 


(Za) • 


(3.2a) 

(3.2b) 

(3.2c) 


By taking partial derivatives of (3.2c) with respect to mi and m 2 , and using the commutation 
of partial derivatives, we obtain the simple result 

^(^ 1 ) _ d{Z2) 

dm 2 dm I ’ 

which can be seen as reflecting some average equilibrium state of the binary under the mutual 
gravitational attraction of its components. This generalizes a similar result established for 
circular motion in Paper I [see Eq. (2.42) there]. 

Next, by taking partial derivatives of Eqs. (3.2), we shall derive relationships between the 
functions M{u,n,ma), L{u,n,ma) and {za){oJ,n,ma), while eliminating R{u),n,ma). First, 
we can express d'^R/dudma in two independent ways by computing the derivative of (3.2a) 
with respect to m^ and the derivative of (3.2c) with respect to cu, yielding 


dL _ d{za) 
dma doj 


(3.4) 


Similarly, we can express d^R/dndma in two ways by computing the derivative of Eq. (3.2b) 
with respect to m^ and the derivative of Eq. (3.2c) with respect to n. Combined with (3.4), 
the resulting expression yields 


dM 

dma 


(Za) 


— U 


djZg) 

doj 


— n 


djZg) 

dn 


(3.5) 


The relationships (3.4) and (3.5) generalize to eccentric orbits Eqs. (4.7) and (4.9) of Paper 
I. These results should prove useful when applied in the context of black hole perturbation 
theory. Indeed, existing gravitational self-force computations already provide numerical data 
for the invariant function {zi){u:, n) for a particle of mass mi on a generic bound (eccentric) 
orbit around a Schwarzschild black hole of mass m 2 3> mi [59, 66, 67]. Combining such 
results with Eqs. (3.4) and (3.5), crucial information regarding the total energy M{u, n) and 
orbital angular momentum L{u,n) of the binary system could be inferred, similarly to what 
has been done in Refs. [15, 18, 20] for circular motion; see Secs. VB and V C below. 

Finally, we can express d‘^R/dudn in two ways by computing the derivative of Eq. (3.2a) 
with respect to n and the derivative of Eq. (3.2b) with respect to cu, yielding 















This expression extends to eccentric orbits the so-called thermodynamic relation commonly 
used in PN theory for quasi-circular orbits (see, e.g., Refs. [54, 68]), or in the construction of 
sequences of quasi-equilibrium initial data for binary black holes and neutrons stars [69-73]. 

Finally, by combining Eqs. (3.2c), (3.4) and (3.5), we can relate the rates of change of 
the radial action with respect to the particles’ masses to the rates of change of the averaged 
redshifts with respect to the radial frequency via 

dR ^ djzg) 
drua dn 


B. First integral relation 

Since Einstein’s equation does not involve any privileged mass scale, the ADM mass M 
must be a homogeneous function of degree 1 in the four variables mi and m 2 , i.e., 

M{XL^R, XR^R, Xmi, Xm 2 ) = X M{L^R, R^R,mi,m 2 ) for any A 7 ^ 0 . Therefore, applying 
Euler’s theorem together with Eqs. (3.1) immediately yields the first integral relation 

M = 2{uL + nR) + ma{za) ■ (3.8) 

a 

This generalizes a similar result established for circular motion in Paper I. Loosely speaking, 
the total mass-energy M is given by the sum of a “azimuthal energy” 2uL, a “radial energy” 
2nR, and the redshifted masses ma{za)- Since Einstein’s equation is nonlinear, it is remark¬ 
able that the total mass-energy of the binary system can be written in such a simple way in 
terms of the individual masses of the particles and other invariant quantities characterizing 
the orbit. 

Alternatively, using the time average introduced in Sec. IIB above, the hrst integral (3.8) 
can be recast in the form 

M = (2{rpr + ipp^)+ '^maZa^ . (3.9) 

a 

However, the different terms in the right-hand side do not possess simple physical interpre¬ 
tations in terms of kinetic energy Ek and gravitational potential energy U. Indeed, in the 
Newtonian limit c~^ —?■ 0, the first term gives 2{rpr -|- g>Pip) = 4(i7k) while the second term 
yields J2a'^a{za) = m — {Ey,) + 2{U). The Newtonian virial theorem implies {U) = —2(Ek), 
and we recover the total energy M = m — (E^) = m + {U)/2 of a Keplerian elliptic orbit. 


IV. VERIFICATION OF THE FIRST LAW 

By making use of available results for the conservative PN dynamics of two point masses, 
we shall now check that the first law of binary mechanics (2.10) and the first integral relation 
(3.8) are indeed satished up to 3PN order. 

We use the result of Ref. [54] for the 3PN-accurate expression of the radial action R{E, L) 
as a function of the reduced binding energy E = (M — m)/p and the dimensionless angular 
momentum L = L/{mp), where m = mi + m 2 is the total mass and p = mim 2 /m is the 
reduced mass. Following Ref. [74], we choose instead to parameterize the orbit in terms of 
the dimensionless variables 

e = -2E, j = -2EL\ (4.1) 


9 




such that in the Newtonian limit e ^ m/a and j ~ 1 — e^, where a and e are the semi-major 
axis and eccentricity of a Keplerian orbit. In terms of the variables e and j, the ADM mass, 
orbital angular momentum and radial action are given by [40, 54, 75] 


1 , ^ 

M = m -e , 

2 ’ 

r 2 [J 
L = m jy \ , 


R = 




35 




15 1 / 

-TV + W + 


15 


-\- — u -\- - 

128 64 128 
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1024 1024 


v + 


15 

1024 


15 

T 

+ 


3 

Tj 

3 

-z/ 


V7 


= + 


35 

T 


-z/ 


1-3/2 


15 15 

- V 

16 16 


V7 


; + 


1024 
231 


-z^"- 


105 109 


125 


- V + 

3 

123 


41 


15 


-TT V H-Z/" 

128 4 


- V H-TT V + 

2 128 


21 


-V 


1-5/2 


(4.2a) 

(4.2b) 


(4.2c) 

1 


p/2 



where u = fi/m = mim 2 /m? is the symmetric mass ratio, such that z^ = 1/4 for equal-mass 
binaries and z/ —)■ 0 in the extreme mass-ratio limit. 

The observables of the orbital motion, namely the periastron-to-periastron period P and 
the angular advance per radial period $ then follow from the first two relations in Eqs. (3.1), 
which can be combined to give P/2ti = dP/dM\L and $/27r = —dP/dL\M- By computing 
the ratio of these expressions we recover uj = ^/P = dM/dL\R. Performing the change of 
variables (M, L) —)■ {e,j) and using the chain rule while computing the partial derivatives of 
(4.2c), we obtain the following 3PN-accurate expansions for the frequencies of the motion: 


n = 


,3/2 


m 


UJ = 


,3/2 


m 


1 + 1 -^ 


+ 


1 + 


+ 


V 

- I £ 


555 15 


11 


T —z^ “t“ 

128 64 128 

9795 1665 


z/ 2 - 


15 


— 3z/ 


^3 


1024 1024 

T05 218 


-V — 


+ 


45 
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- V 


105 
1024 
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- T\ V H- V" 

64 2 


15 


z/^ + 


1-3/2 


255 135 


- -^z/ + 


15 


-V 


Vj 


+ o{e^) 


(4.3a) 
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'2685 225 


105 2 
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(4.3b) 
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Equation (4.3a) agrees with Eq. (7.7a) of Ref. [74], while Eq. (4.3b) can easily be recovered by 
combining their Eqs. (7.7a) and (7.7b). At Newtonian order, we have u = n, or equivalently 
A$ = 0, which is to say no periastron advance. As is well known, bound Keplerian orbits 
are closed ellipses. 

The formulas (4.2c) and (4.3) were derived from the expression for the 3PN binary Hamil¬ 
tonian, in the center-of-mass frame, in ADM-TT coordinates [54]. The coordinate-invariant 
relations {za){e,j) were, on the other hand, recently computed, up to 3PN order, by following 
an entirely different route. Indeed, the near-zone metric gapiyi) = dagit, Yi) evaluated at the 
coordinate location yi{t) of the particle 1 was computed up to 3PN order, in harmonic coor¬ 
dinates [43]. Akcay et al. [ 66 ] then obtained the redshift zi = {—gapiyi) vfviY^'^ of this par¬ 
ticle by contracting the near-zone metric with the coordinate velocity vf = dyf /dt = ( 1 , v\). 
Specializing the resulting expression to the center-of-mass frame, using the 3PN generalized 
quasi-Keplerian representation of the motion [76], and performing an average over one radial 
period, they obtained the following 3PN-accurate expression for the average redshift:® 
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(4.4) 


Here, A = (m 2 —mi)/m = (1 — 4 z /)^/2 jg reduced mass difference. (We assume mi < m 2 .) 
The expression for {z 2 ){e,j) is easily found by setting A —)■ —A in Eq. (4.4). 


5 Reference [66] computed instead the average oiu\ = 1/zi with respect to proper time (denoted (U) there); 
we simply have ( 2 : 1 ) = 1/(C/). 
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Now it is straightforward to check that the 3PN results (4.2)-(4.4) do indeed obey the first 
integral relation (3.8), as well as the partial differential equations (3.2) with the substitutions 
{d/du,d/dn) —> {d/de,d/dj). This verification provides a powerful check of the intricate 
calculations that resulted in the expression for the 3PN Hamiltonian in ADM-TT coordinates 
[54], of the 3PN near-zone metric in harmonic coordinates [43], and of the hrst law itselfy. 
Moreover, if the redshift {zi) were computed up to 4PN order, the verihcation of Eq. (2.10) 
would provide a powerful test of the recently derived 4PN binary Hamiltonian [35-37]. 

Finally, we discuss the reduction to the case of circular orbits; degenerate orbits for which 
the frequencies uj and n are not independent. By dehnition, a circular orbit has a vanishing 
radial action integral. Setting i? = 0 in Eq. (4.2c) yields a relationship between e and j (or 
equivalently between M and L) that reads, up to 3PN order [54, 74], 


j = 1 + 


9 V 
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e + 


81 1/2 

-2z/ H- 

16 16 


/2835 
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192 


7699 41 , 

-TT 


392 32 


V 




1 / H-^-U + oie'^). 

2 64 / ^ ^ 


(4.5) 


Substituting for this j into (4.3b), and inverting the resulting formula, we recover the well- 
known 3PN-acurate expression for the total mass-energy M as a function of the circular-orbit 
frequency oj. By inserting this expansion into Eqs. (4.5) and (4.4), we also recover the known 
3PN results for L{uj) and .^ 1 ( 0 ;) for circular orbits; see, e.g., Eqs. (2.35)-(2.37) of Ref. [15]. 


V. APPLICATIONS OF THE FIRST LAW 

In this section we discuss several applications of the first law of binary mechanics, starting 
with the gravitational-wave driven adiabatic inspiral of compact binary systems (Sec. V A). 
We then explain how this relation could be used, together with existing and/or forthcoming 
GSF calculations (Sec. VB), to get strong-field information about the conservative dynamics 
of compact binaries, especially regarding the binding energy, angular momentum and radial 
action of such systems (Secs. V C and VD). Finally, we illustrate how our results can be used 
to inform the conservative part of the EOB dynamics beyond circular motion (Sec. VE). 


A. Gravitational-wave driven adiabatic evolution 

The variational first law relates the total energy, orbital angular momentum, radial action 
integral, and averaged redshifts of two physically distinct binary systems under small changes 
of the orbital frequencies cu, n and of the particles’ masses mi, m 2 . In the problem of building 
template waveforms for inspiralling compact-object binaries, we wish to follow the evolution 
of a single system as it gradually inspirals under the effect of gravitational radiation reaction. 
In this case, the quantities M, L, R and tUq are no longer constants of the motion. 

Nevertheless, if the characteristic timescale Ti-.r. of gravitational radiation reaction is much 
larger than the typical orbital timescale P, then the time evolution of M{t), L{t), R{t) and 
ma(t) is well approximated as an adiabatic process. We may then identify the two physically 
distinct systems compared in the first law with two nearby states of a single binary system. 
This approximation is commonly adopted while computing sequences of quasi-equilibrium 
initial data for binary systems of black holes and neutrons stars [69-73]. The variations 5M 
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and 6L appearing in Eq. (2.10) are then interpreted as the secular changes in the mechanical 
energy and angular momentum of the binary during an interval P < -C Tr r,, and similarly 
for 6R and Srria- Dividing (2.10) by 6t and averaging over a radial period, the average rates 
of change of M(t), L(t), R{t) and ma(t) must then obey 

{M) = uj (L)+n{R) {ma) . (5.1) 


Notice that the frequencies uj, n and {za) were factored out from the averaging because they 
are already averaged over one radial period and they do not vary over a dynamical timescale. 

We may then appeal to the usual argument of balance of energy and angular momentum, 
(M) = —P and (L) = —Q, where the gravitational-wave fluxes of energy P and angular mo¬ 
mentum Q can be computed from the far-zone gravitational held [77, 78]. These (heuristic) 
balance equations are motivated by the (exact) Bondi-Sachs mass-loss formula [79, 80]; see 
e.g. Sec. II E of Paper I for a discussion. Moreover for a neutron star the conservation of the 
baryonic mass implies (rha) = 0, while for a black hole tidal heating [81, 82] is responsible 
for an increase in the irreducible mass, the average rate of which is given by the hux of 
energy Ra through the horizon: (rha) = Ra [83, 84]. Therefore, Eq. (5.1) becomes 

n{R) = -R + ujQ - '^ea{za)Ra, (5.2) 


where £« = 0 for a neutron star and £« = 1 for a black hole. 

For a nonspinning black hole, the tidally induced hux of energy Ra through the horizon 
is a small, 4PN ehect relative to the leading-order (Newtonian) huxes P and Q [85], and will 
thus be neglected here. Using the known expressions for the leading-order huxes of energy 
and angular momentum [86, 87], we get for any binary system of compact objects 


(R) = —— 
5 


(1 - e2)"/" 


73 


37 


1 + wre + —e 


24 


96 




where x = (mcu)^/^ is the usual, dimensionless, frequency-related PN parameter, and e is the 
eccentricity parametrizing the Keplerian orbit. At this order of approximation, j = 1 — 
and n = oj. Now let /(e) denote the factor in curly brackets in Eq. (5.3). It is positive for all 
0 < e < 1 and has the asymptotic behaviors /(e) ~ y e^ and /(e) ~ (^ “ e)“^/^ in the 

limits where e —?■ 0 and e ^ 1, respectively. Hence (P) < 0 and we recover the known result 
that, in the weak-held regime, gravitational radiation reaction decreases the noncircularity 
of the orbit, all the more so if the eccentricity is large [87]. Of course (P) vanishes in the limit 
e —)■ 0. For circular orbits, the “adiabatic hrst law” (5.2) shows that the huxes of energy 
and angular momentum must be proportional: R = ujQ. By using available PN expressions 
for these huxes, this relationship can be checked up to 3PN order [74, 78]. 

For completeness, we note that the average rates of change (h) and (c5) of the frequencies 
entering Eq. (5.1) are given by 


(n) 

(D) 


dn 

du 


^ dn ^ 


dL 

_ doj V—\ 


dn 

drua 

doo 

dnia 


Pa, 

(5.4a) 

Pa, 

(5.4b) 


where the partial derivatives with respect to M, L and can easily be computed by means 
of the chain rule, together with Eqs. (4.1) and (4.3). 
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B. Conservative dynamics beyond the geodesic approximation 


The derivations of the hrst law of binary mechanics given in section II did not rely on any 
small-velocity or weak-held expansion. However, they both relied on several key assumptions 
that are typically met by PN spacetimes, such as (i) a well-dehned conservative/dissipative 
split of the binary dynamics, (ii) the existence of an autonomous Hamiltonian describing the 
conservative dynamics, (iii) a point-particle description of the compact objects, and (iv) the 
implicit use of a regularization scheme to subtract off the divergent self-helds of the particles. 
Therefore the applicability of the formula (2.10) to modelling compact-object binaries in the 
strong-held regime, for which the PN approximation breaks down, is not guaranteed. 

However, several results suggest that Eq. (2.10) does indeed hold in this context. First, 
by making use of the circular-orbit hrst law of Paper I, the authors of Refs. [18, 19] could 
recover the exact (numerical) value of the shift of the Schwarzschild ISCO frequency induced 
by the conservative part of the GSF [88, 89]. Moreover, Eq. (2.10) is formally identical to 
the (nonspinning limit of the) hrst law-type relations of Refs. [90, 91], which were derived 
in the context of black hole perturbation theory. In particular, Isoyama et al. [91] have 
devised a Hamiltonian formulation of the dynamics of a self-gravitating particle subject 
to the conservative GSF, for bound orbits in a Kerr background; the restriction of their 
hrst law relationship to circular equatorial orbits was used to compute the GSF-induced 
frequency shift of the Kerr ISGO [92]. Independently, Vines and Flanagan [93] have recently 
proved that, for generic stable bound orbits in a Schwarzschild background, the dynamics of 
a pointlike object subject to the conservative piece of the osculating-geodesic-sourced GSF 
[94] is Hamiltonian and integrable. 

Thereafter, we shall thus assume that the hrst law (2.10) (with 6 m 2 —)■ 0) can be applied, 
together with perturbative GSF calculations, to obtain new information about the conserva¬ 
tive dynamics of nonspinning compact binaries, even in a strong-held regime. Alternatively, 
one might adopt the viewpoint that the formula (2.10) defines, in the perturbative context, 
some physically motivated notions of binding energy, angular momentum and radial action. 


C. Binding energy, angular momentum and radial action 

Before we proceed to discuss these applications, let us hrst establish a few useful formulas. 
Introducing the variable Ai = M — uL — nR, the relationships (3.4)-(3.7) can be combined 
to express M, L, R and {za) solely in terms of Ai and its partial derivatives with respect to 
uj, n and rua as 


AI = 

., dAi dAi 

ouj on 

(5.5a) 

T 

dAi 

(5.5b) 


du ’ 

R = 

dAi 

(5.5c) 

dn 

Za) = 

dAi 
drua ’ 

(5.5d) 


Thereafter it will prove convenient to introduce the dimensionless quantities E = [M—m) //i, 
L = L/{m^), R = R/{mfi) and A4 = {Ai — m)//i, and to perform the change of variables 
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{mi,m 2 ,n,uj) —> {m,i>,n,uj), where h = mn and Cj = moj are dimensionless “versions” of 
the fundamental frequencies. In particular, since for dimensional reasons Ai{n,uj, u) cannot 
depend on the total mass m, Eq. (5.5d) can be recast in the form 


{Za) = 1 + - (1 ± A - 4i/) 


M + v 


dM 

dv 


^dM .dM 

M + oj + n 


dCj 


dn 


(5.6) 


where the plus sign (respectively, minus sign) stands for particle 1 (respectively, particle 2), 
and we recall that A = (m 2 — mi)/m = (1 — is the reduced mass difference. 

We now restrict to the extreme mass-ratio limit q = mi/m 2 1, and denote by (z) = {zi) 
the average redshift of the small body. Although it is common, in the context of perturbation 
theory, to use the small mass ratio q as an expansion parameter. Refs. [14, 18, 62, 64, 95-100] 
showed that it is more advantageous to use the symmetric mass ratio v = q/ (l-|-g)^ instead.® 
(Formally, q = u + Hence, to leading order beyond the test-particle approximation, 

we consider the expansions 


E{o) + V E{i) -b C>(z/^), 

(5.7a) 

L(o) -b -b Oiy ^), 

(5.7b) 

R(o) -b , 

(5.7c) 

A4(o) -b nM(i) + 0{u ‘^), 

(5.7d) 

(z)(o) -b n {z)(i) + 0{i /), 

(5.7e) 


where all the quantities are functions of the dimensionless frequencies n and A. The variables 
with a subscript (0) correspond to the geodesic values, while those with a subscript (1) denote 
the conservative GSF corrections (at hxed frequencies).'^ 

At zeroth order in the (symmetric) mass ratio z/, Eq. (5.6) implies (2()(o) = 1 + A4(o), such 
that Eqs. (5.5a)-(5.5c) yield the following relations between the leading-order contributions 
to the average redshift, binding energy, angular momentum and radial action integral: 


E{o) 

L(o) 


(^)( 0 ) 


— u 


du 


duj 


R(o) 


— n 


d{z)(o) 

dn 


- 1 , 


dn 


(5.8a) 

(5.8b) 


Adjusting notations, these formulas agree with Eqs. (A2)-(A3) and (A7) of Ref. [66]. At the 
next-to-leading order, (5.6) implies (z}(i) = 2M(i)—2M(o)(dM(o)/dcj)+ n (dM(o)/dn), 
from which we deduce that the GSF contributions to the binding energy, angular momentum 
and radial action integral are given in terms of the GSF contribution to the average redshift 
(as well as geodesic quantities) by 


^( 1 ) = 


/ / \ -- ^W(i) 


— n 


dn 


+ 2F(o) 


+ A 


2 ' 9 ^(^)(o) , . d‘^{z){o) , .2 

-f- 2u:n + n 


dCj"^ 


dCjdn 


d‘^{R)(o) 

dn^ 


(5.9a) 


® Note that this parameter appears naturally in the PN expansions of quantities that are symmetric under 
the exchange 1 •<->■ 2 of the bodies’ labels; see e.g. Eqs. (4.2c) and (4.3). 

^ In particular, l-l-iil(Q) and E(o) coincide with the usual conserved (specific) energy and angular momentum 
associated with the stationarity and axisymmetry of the background Schwarzschild geometry. 
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L(i) 

R{1) 


1 f d{z)(i) d{z)(0) _ ^ d‘^{z)io) _ ^ d‘^{z)(0) \ 

2 \ 5a; dCj ^ dhdoj ) ’ 

1 ( d{z)(i) d{z)io) _ ^ 5^(^)(o) _ . 5^(^)(o) \ 

2\dh dh ^ dfi^ ^ dudh ) 


(5.9b) 


(5.9c) 


These relationships generalize Eqs. (4) of Ref. [18] to generic bound (eccentric) orbits. Hence, 
if the function {z)(i)(a;,h) is known, the leading-order corrections £'(i)(a;,h), L(i)(a;,h) and 
R(i)(a;,h) to the geodesic values can be computed for any bound orbit. Such GSF data for 
( 2 ;)(i) is already available for a range of orbits with eccentricities 0 < e < 0.4 [59, 66, 67]. 

The relationships (5.9) could easily be generalized up to second (or higher) order, yielding 
the expressions for the 0 {v‘^) contributions to the binding energy, angular momentum and 
radial action integral in (5.7), say 77(2)(h, a;), L{ 2 ){n,Cj) and i?( 2 )(h,a;), in terms of the second 
order GSF contribution to the average redshift, say ( 2 ;)( 2 )(h,a;), and its frequency derivatives 
(as well as geodesic and hrst-order GSF quantities). Thanks to recent progress in formulating 
the second-order GSF [101-107], numerical results for the second-order contribution Z( 2 )((^) 
to the redshift for circular orbits should soon become available [108]. Gombining these data 
with the hrst law of binary mechanics would then yield the fully relativistic expressions for 
the functions 77(2) (a;) and L( 2 )(a;), to be added to the geodesic and hrst-order contributions 
that were obtained in Ref. [18]. 


D. Schwarzschild separatrix and singular curve 

Interestingly, Fqs. (5.8) and (5.9) could be used to explore some conservative GSF effects 
on the motion of a self-gravitating particle on a bound orbit around a nonspinning black hole. 
For instance, one could compute the GSF-induced shift in the location of the Schwarzschild 
separatrix, the curve that separates—in the relevant parameter space—between stable bound 
(eccentric) orbits and unstable (plunging) ones [109]. (Along the separatrix, the radial period 
P and the angular advance per orbit <I> both diverge, while their ratio $/P remains hnite; 
therefore these orbits are characterized by a vanishing radial frequency n at a hxed a;.) Since 
all of the orbits that lie along the separatrix are marginally unstable, this calculation would 
generalize to eccentric orbits that of the shift of the Schwarzschild ISGO frequency induced 
by the conservative GSF [18, 19, 88, 89, 92]. Especially interesting would be the calculation 
of the GSF-induced shift in the frequency of the Schwarzschild innermost bound stable orbit 
(IBSO). Indeed, the ISGO and the IBSO are the end points of the separatrix. 

However we note that, in order to compute the binding energy and angular momentum 
of such marginally unstable orbits, it may not be necessary to rely on eccentric-orbit GSF 
data for (z)(i)(a;, h), in the limit where h —)■ 0 at hxed u). Indeed, in the test-mass approxi¬ 
mation, it is known that to each orbit along the separatrix corresponds an unstable circular 
timelike geodesic (homoclinic orbit [110]) that shares the same binding energy and angular 
momentum. Assuming that this property still holds while taking into account the effect of 
self-interaction, one could simply combine the (circular-orbit restriction of the) hrst law of 
binary mechanics with the circular-orbit GSF data for Z(^i){uj) that is provided in Table IX 
of Ref. [19]. 

Another key property of bound timelike geodesic orbits around a Schwarzschild black hole 
is the existence of a singular curve in the parameter space, along which the Jacobian matrix 
of the transformation (77(o), 7^(0)) •H- (a;, n) becomes singular [59]. This strong-held feature is 
closely related to the existence of a separatrix, and is responsible for the recently discovered 
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phenomenon of isofrequency pairing [111], namely, the existence of physically distinct orbits 
(having different i?(o) and -h(o)) that share the same frequencies u and n. Going beyond the 
test-particle approximation, we hnd using Eqs. (5.7)-(5.9) that the Jacobian determinant of 
the transformation {E, L) -H- (a), n) is given by 


d{E,L) _ . f f d^z) (0) \ ^ {z) (p) {z) (p) \ 

d{u,n) ^ V dudn J dn?' j 


dudn \ dudn ^ du'^dn ^ dudn^ ) 


2 du"^ \ dn^ dudfi^ dh^ ) 


2 dn'^ \ du"^ du^ du'^dn J 


+ 0{u^). (5.10) 


In the test-mass limit —?■ 0, the Schwarzschild singular curve u = Usin) along which J —)■ oo 
can be computed from the condition [{d‘^{z)(^o)/dudn)‘^ — {d‘^{z)(^o)/du‘^){d‘^{z)^o)/dn‘^)] 0. 

Moreover, once the second-order partial derivatives of {z)(^i){n,u) evaluated along the curve 
u = Us{n) are known, it will become straightforward, thanks to Eq. (5.10), to compute the 
GSF-induced shift in the location of the Schwarzschild singular curve. 

Carrying out the applications outlined above requires obtaining GSF data for (z)(i)(h, a)) 
(as well as its partial derivatives) in a region of the parameter space that remains challenging 
for state-of-the-art GSF codes, namely for strong-held geodesic orbits that lie deep in the 
“zoom-whirl” [112] regime.® Nevertheless recent work suggests that such computations could 
become tractable in the foreseeable future, for instance by using a method relying on Green 
functions and worldline integrations [114]. 


E. Noncircular conservative EOB dynamics 

As pointed out in Ref. [32], GSF results can be used to inform the EOB model, including 
for strong-held orbits, by providing the exact expressions for the 0 {u) contributions to the 
metric potentials that enter the conservative piece of the EOB dynamics. By construction, 
the conservative part of the EOB dynamics derives from the Hamiltonian [11] 

-^EOB = m sjl + 2v — l) , (5-11) 

where m = mi -|-m2 is the total mass of the binary, while hfes = TJes/h denotes the “ehective 
Hamiltonian” of an “ehective particle” of mass /i = mim^/m that follows a timelike geodesic 
(modulo quartic and higher-order terms in the momentum) in the “ehective metric” 

dSeff = g%{,x) da;“da;^ = —A{r] v) df^ J- i?(r; v) dr^ J- (d6*^ J- sin^ 9 . (5.12) 


Unbound orbits close to the separatrix also play an important role in a scenario that explores the possibility 
of overspinning a nearly extremal Kerr black hole by means of a particle plunging from infinity [113]. 
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This effective metric is a deformation of the Schwarzschild geometry of a black hole of mass 
m, with deformation parameter v = fi/m, such that A = B~^ = 1 — 2 m/r in the test-mass 
limit —?■ 0. The effective Hamiltonian entering the EOB Hamiltonian (5.11) reads [115, 116] 


HefT(r,pr,Lo) 



r2 


pI 




+ Q{r, pr] y) 


(5.13) 


where Lq is the conserved EOB angular momentum and Q{r,Pr] v) is a function that controls 
the deviation from geodesic motion in the effective metric Up to 3PN order, we have 

Q = (8 — 6z/) vp\lr^ where, for convenience, we introduced the rescaled variables f = r/m, 
Pr = Pr!p and Q = Q/p^- The 4PN contribution to the EOB potential Q, and those entering 
the expressions for A and B, have recently been computed in Ref. [116]. Following Ref. [32], 
we shall assume that the fact that Q depends solely on pr and r, and vanishes (at least) like 
Pr when Pr —)■ 0, remains true at higher PN orders; see also [20]. Moreover, we shall restrict 
our analysis to mildly eccentric orbits, for which the contribution 0 {p^) (and higher orders) 
to Q can be neglected. 

Now, in the extreme mass-ratio limit z/ 1, the potentials A{u; z/), B{u; u) and Q{u,Pr] v) 
that enter the conservative EOB dynamics can be expanded as 

A = 1 — 2u + u a{u) + 0{u ‘^), (5.14a) 

D = l + iyd{u) + 0{u^), (5.14b) 

Q = z/ q{u) Pr + , (5.14c) 


where m = 1/f is a dimensionless measure of the gravitational potential, and D = {AB)~^. 
For circular orbits, all of the information about the conservative EOB dynamics is encoded in 
the metric coefficient = —A. By combining the hrst law of Paper I with GSF data for the 
redshift variable. Refs. [18, 20] could compute the function a{u) for all 0 < m < 1/5. Making 
use of additional GSF data for circular orbits, all the way down to the Schwarzschild light¬ 
ring, Akcay et al. [19] then completed the determination of the exact linear-in-z/ contribution 
to the EOB potential A{u; u) for all 0 < m < 1/3. Moreover, Damour [32] showed that, for 
slightly noncircular orbits, the GSF contribution to the invariant relationship K{oj) is related 
to a linear combination of the functions a{u) and d{u). This GSF effect was later computed 
numerically in Ref. [60], which allowed to determine the linear-in-z/ contribution to the EOB 
potential D{u; u) [and hence B{u; z/)] for all 0 < tz < 1/6 [19, 20]. However, the function q{u) 
that encodes the C>(z/p^) contribution to the EOB potential Q{u,pr; z/) has, so far, remained 
entirely unconstrained (except for the leading 3PN and subleading 4PN terms when tz —?■ 0). 
We shall now show how, by combining Eqs. (5.9) with GSF data for the relation (z)(i)(fz, a)) 
(for mildly eccentric orbits), the function q{u) can be determined. 

Having neglected terms 0{p^) (and higher) in the EOB potential Q, and working to linear 
order in z/, the effective EOB Hamiltonian (5.13) squared can be expressed as a quadratic 
polynomial in the radial momentum squared as 


Pr, Lq) = A{u] Z/) (1 -F A{u] 1>)D{U] i')pI + u q{u) p. 


(5.15) 


Inserting the expansions (5.14), and solving for the radial momentum yields Pr as a function 
of tz = 1/f and the conserved quantities Sq = Hgs and Lq = LQ/{mp). Splitting the resulting 
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expression into 0{iy^) and 0{iy) contributions, i.e., pr = we find 


p‘;^\u;io,Lo) 

2p^^\u;So,Lo) 


£^-{l-2u){l + u^Ll) 


P2 

do) n - 


1 -2n 
a{u) 


(1 - 2n)^ 


( d{u) + 


a[u) 

l-2u 


+ m‘‘ 


q{u) 

l-2u 


(5.16a) 

(5.16b) 


As expected, Eq. (5.16a) coincides with the well-known expression for the radial momentum 
of a test particle of energy Eq and angular momentum Lq in orbit around a nonspinning black 
hole of mass m. The correction term (5.16b) involves, additionally, the known functions a{u) 
and d{u), as well as the unknown function q{u) that we intend to constrain. 

For bound orbits, the turning points of the radial motion correspond to the two smallest 
real, postive and finite roots of the equation Pr{u] Eq, Lq) = 0, say = l/r.,. and u+ = l/f_. 
To linear order in z/, these can be written as 



(5.17) 


where the geodesic values are known, for any given Eo,Lo, by solving pf ''*; Eg, To) = 0. 
Expanding pr(u; Eg, Lq) = 0 to linear order in z/, the mass-ratio corrections are found 

to be solutions of u^\dpi^^/du)(u^^) +pi^\u^^) = 0, the expression of which requires the 
knowledge of the unknown values for g(u±^). However, as we shall see shortly the knowledge 
of will not be necessary to constrain the function g(u). 

For any given Eg and Lg, we may now integrate the radial momentum (5.16) over a radial 
period to obtain the reduced EOB radial action (recall that u = 1/f) 


^EOB = — (p Prdf = 



(5.18) 


Using the formula (5.17), this integral splits into three contributions, the ranges of integra¬ 
tion of which vary from to from to u^\ and from u® to 

Expanding to linear order in z/, the first and last integrals will contribute terms proportional 
to z/M^^pr°^(M^^) and z/M^^pr°^(M+^), respectively. But these vanish by virtue of the defi¬ 
nition of the radial turning points in the test-mass limit, such that the EOB radial action 
(5.18) can be computed by integrating over the geodesic interval between u^^\Eg, Lg) and 
m+^(To, Lg). Hence, for any given Eg and Lg, the reduced EOB radial action may be written 
in the form i?EOB = R(o) + ^ where i?(o) is the radial action for a test particle 

orbiting a Schwarzschild black hole of mass m, and 

1 /■“+' - du 

AR = - pld(u-, Eg, Lg) — . (5.19) 

TT y„(o) 


Next, recall that the EOB model is (by construction) built upon the identification of the 
“on-shell” Hamiltonian (5.11) with the total mass-energy of the binary, TfEOB = M, as well 
as the identification of the EOB angular momentum entering the effective EOB Hamiltonian 
(5.13) with that of the real two-body system, i.e. Lg = L [11]. Therefore, in terms of reduced 
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variables, we find by inverting (5.11) with M = m + fiE that the EOB and real radial actions 
must obey 


Reob{^o, Lq) = R{E, L ), where 


£^ = l + E + (z//2) E2 ^ 
Lo = L. 


(5.20) 


The perturbative equalities R = REOB + {i^/‘^)E^{dR/dE) + 0{u'^) and .Reob = R{o) + J^ ^R+ 
0{u‘^) can be combined, together with the relationship OR/dE = n~^ and OR/dL = —<|)/27r, 
as well as Eqs. (5.7a) and (5.7b) to expand R{E,L) at fixed frequencies n and ca, to obtain 

{ 1 E'^ 1 

x(£,i| - + ^ + Afl + 0(v^). (5.21) 


Combining Eqs. (5.8a) and (5.9) yields i7(i) — = hi?(i) + |((^)(i) + (^)(o) + -^(o) ~ l)- 

Hence, using Eq. (5.7c) we note that the sum of the first two terms in the right-hand side of 
Eq. (5.21) cancels out the left-hand side. Since the 0{i>) contribution must vanish identically, 
we obtain the key identity 

2n AR + (^)(i) + (^)(o) + E{o) + Efo) = 1 • (5.22) 

In summary, we have shown that—for any given fi and Cj —the definite radial integral (5.19), 
which involves the functions a{u), d{u) and q{u), is equal to a known function of the funda¬ 
mental frequencies (assuming that the GSF contribution {z)(i){fi, u) is known). In particular, 
Eq. (5.22) shows that, for mildly eccentric orbits, all of the information about the linear-in-z/ 
contributions to the EOB potentials is encoded in the function (z)(i)(h, cD). Using Eqs. (5.16) 
and (5.22), we can now easily deduce explicit expressions for a{u), d{u) and q{u). 

To do so, it is quite convenient to parameterize the orbit in terms of a “semi-latus rectum” 
p and an “eccentricity” e, instead of the frequencies h and u, and to perform an expansion in 
the limit e —?■ 0, which is well adapted to the mildly eccentric orbits that we are considering. 
Following Darwin [117], we parameterize the radial motion using the “relativistic anomaly” 
X via 

P 

r{x) = w-- u{x) = n (1 -t- e cos x), (5.23) 

1 + e cos X 

where v = 1/p is the inverse semi-latus rectum and 0 < e < 1 the eccentricity. In terms of the 
parameters {v, e), the boundaries of the integral (5.19) simply read u® = v (l±e). Adjusting 
notations, the expressions for h{v,e), oj{v,e), {z)(o){v,e), i7(o)(n,e) and L(o)(n,e) are given, 
e.g., in Eqs. (2.4)-(2.10) of Ref. [66]. All of these relationships can be computed analytically, 
as perturbative expansions in powers of e. Moreover, performing the change of variable u 
X and substituting for (5.23) in Eq. (5.19) [with (5.16)], the integral AR(i7(o)(n, e), T(o)('y, e)) 
can also be computed analytically, in the small-eccentricity limit. The result, which we have 
obtained up to is too cumbersome to be displayed here. Finally, in the small-e limit, 

the GSF contribution to the generalized redshift can be expanded as 

(2)(i)(n,e) = z^i){v) + {z)f^){v) + ^ (^)p)(n) + o(e^), (5.24) 

where we used the notations Z(i)(n) = hme^.o {z){i){v^e) , (-s)(i)(n) = lime^.o 9^(-s)(i)(n, e)/5e^, 
and {zY^^^{v) = \im.e.^Qd‘^{z){i){v^e)/de^. The terms 0{e) and 0{eY must vanish, as can be 
checked explicitly from the absence of such contributions in the identity (5.22). 
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Equations (5.16) show that ~ e and ~ 1/e in the limit where e ^ 0. Therefore, 
the leading-order contribution in Eq. (5.19) is C>(e°). Similarly, the variables h, (2;)(o), (^)(i) 
and E(o) all contribute at leading O(e^). At that order, the identity (5.22) implies 


a{v) = \/l — 3n Z(^i){v) — v 



1 — 4n \ 


(5.25) 


This formula agrees with the result obtained in [20] (see (2.14) therein). While the derivation 
given in Ref. [20] required integrating an ordinary differential equation, our derivation yields 
in a straightforward manner the algebraic relation between a and Z{iy 

At the next O(e^), the identity (5.22) provides a relationship between the functions d(n), 
a(n), a'{v) = da/dn, a''{v) = d^a/dn^ and (2)(i)(n). Substituting for Eq. (5.25), we hnd 


d{v) = 


1 - 


(7 - + 45n2) 

2 (1 - 3n)®/^ (1 - 3n)^/^ 


Ri)(^ 


2 - f n + 101^2 


132^3 


(1 — 6v) \/l — 3v 


hi) 


t) 


n (1 - 2v) y/1 - 3v (1 - 2v) \/l - 3v , . 

-%)W+ „(i-6„) 


2(1 -6n) 


(5.26) 


By making use of Eqs. (4.49)-(4.50) and (Bl) of Ref. [66], one can determine the expressions 
for . 2 ( 1 ) (n) and (.2)(i)(n) in the weak-held limit n 1, and then check that Eq. (5.26) provides 
the correct 3PN expansion for the linear-in-i^ contribution to the EOB potential D, namely 
d{v) = 6n^ -|- 52n^ -I- Oipj'^). Moreover, a comparison of the formula (5.26) with Eqs. (5.21)- 
(5.25) of Ref. [32] (together with (5.25) here) yields a relation between the GSF contribution 
to the ratio W = {n/uY (denoted p in Ref. [32]), and the quantities Z(i) and {zYf^-^y It would 
be interesting to use GSF data for p, Z(^i) and to check this prediction. 

At the next 0{eY, we obtain a relation between q{v), d{v), d'{v), d"{v), a{v), a'{v), a"{v), 
a"'(n), a""(n) and (z)®^^(n). Substituting for Eqs. (5.25) and (5.26), we get 




7 (1 - 2vY 99 

- QvY\/1 — 3n V 4 


(1 


6n(l - 6n)’Vl - 3n V ^ 

)vYl-Qvf ^ ^ n(l-3n)' v 

372085 o 467057 . 185935 , 243789 ^ 269793 . A , 

—-—n -I-64188n I z)] 
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- V H- 

3 48 
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T 7a 7= 
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(5.27) 








































We checked that this expression reproduces the known 3PN result for the linear-in-z/ contri¬ 
bution to the EOB potential Q, namely q{v) = -|- 0{v^). Using GSF data for zp) and p, 

the authors of Refs. [19, 60] devised accurate global hts for the functions a{v) and d{v), that 
are based on simple analytic models. While a{v) is known for all 0 < u < 1/3, the function 
d{v) has only been determined over the range 0 < u < 1/6. Using additional GSF data for 
( 2 ;)®]^) and together with Eq. (5.27), one could similarly construct a global £t for q{v). 

Finally, we note that the knowledge of the GSF contributions to the invariant relationships 
M{n,u)) and L{n,u) [from Eq. (5.9) and GSF data for ( 2 ;)(i)(h,a;)] could be used to determine 
the exact, linear-in-z/ contribution to the EOB potential Q{u,pr] v) for generic bound orbits, 
and not merely the function q{u) that encodes the 0{p^) contribution therein [32]. This task 
is left to future work. 
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